Self-vacancies in Gallium Arsenide: an ab initio calculation 



O 

o 

(N 
D 



C/2 



Fedwa El-MellouhQ and Normand Mousseai:OI 

Departement de physique and Regroupement quebecois sur les materiaux de pointe, 
Umversite de Montreal, CP. 6128, suae. Centre-ville, Montreal (Quebec) H3C 3J7, Canada 

(Dated: February 2, 2008) 

We report here a reexamination of the static properties of vacancies in GaAs by means of first- 
principles density-functional calculations using localized basis sets. Our calculated formation en- 
ergies yields results that are in good agreement with recent experimental and ab-initio calculation 
and provide a complete description of the relaxation geometry and energetic for various charge state 
of vacancies from both sublattices. Gallium vacancies are stable in the 0, — , —2, —3 charge state, 
but Vq^ remains the dominant charge state for intrinsic and n-type GaAs, confirming results from 
positron annihilation. Interestingly, Arsenic vacancies show two successive negative-U transitions 
making only +1 , —1 and —3 charge states stable, while the intermediate defects are metastable. 
The second transition ( — / — 3) brings a resonant bond relaxation for V^^ similar to the one identified 
for silicon and GaAs divacancies. 

PACS numbers: 61.72.Ji, 71.15.Mb,71.15.Pd, 
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I. INTRODUCTION 

Native point defects are involved in virtually every pro- 
cess during which an atom incorporated in the lattice 
of a semiconductor migrates toward another lattice site. 
This diffusion mediated by point defects is responsible 
for a number of important effects, for instance, those en- 
countered during fabrication of microelectronic devices. 
It is not surprising therefore that point defects in semi- 
conductors have been extensively studied using ab initio 
techniques. Some work still remains to be done, however, 
especially in the case of alloyed semiconductors. 

If Ga vacancies in GaAs are relatively well understood 
— because most dopants used in technology (Si donor, 
Zn, Be, and Mg acceptors) occupy the Ga sites — much 
less is known about the As vacancy. The introduction of 
carbon as a prospective As-site acceptor causes a grow- 
ing interest in this defect, however. There is therefore an 
abundant literature on calculations of formation energies 
of point defects in GaAs which is nicely reviewed and 
summarized in a paper of Deepak et alX. Since forma- 
tion energies are difficult to measure, calculations, are the 
primary method for obtaining these values. However be- 
cause of the assumptions and approximations taken into 
account, the reported values in the literature during the 
last two decades differ greatly from paper to paper. 

Most calculations, were done for small supercells. 
These calculations take advantage of error cancellations 
to obtain energy differences that are more precise than 
the total energies themselves. Here we repeat these cal- 
culations for all possible charge states for both Ga and 
As vacancies using the strict convergence criteria and a 
large simulation cell. 

The outline of the paper is as follows. We first de- 
scribe the method for defect calculation in section ^1 
then we turn to the convergence tests made to set up 
the methodology. We focus on the effect of the fc-point 
sampling, system-size and completeness of the basis on 
the formation energy and structural properties of the de- 



fects. Section IIIII deals with these effects on the most 
important charge states like V^^ and Vq^ . In section HVI 
we present and discuss the most converged results using 
localized basis sets compared to previous results obtained 
from theory and experiment. 



II. SIMULATION DETAILS AND 
METHODOLOGY 

A. Total energy calculations 

The total energies for this work are evaluated using 
SIESTAiSi^, a self-consistent density functional method 
(DPT) within local-density approximation (LDA). Core 
electron are represented by the standard norm-conserving 
TrouUier-Martins pseudo-potentials^ factorized in the 
Kleiman-Baylander form^ and the one-particle problem 
is solved using linear combination of pseudo-atomic or- 
bitals (PAO) basis set of finite range. These orbitals are 
strictly localized and represent well the local electronic 
densities; few of them are therefore needed, decreasing 
considerably the computational costs by comparison with 
standard plane-waves calculations. The main drawback 
of this approach, however, is the lack of a systematic 
procedure to ensure a rapid variational convergence with 
respect to the number of basis orbitals and to the range 
and shape of each orbital. Consequently, while extending 
plane-wave basis sets is trivial, some efforts are needed 
to prepare unbiased pseudo-atomic basis sets (see, for 
example, Ref||6ij3). 

In this work, we use the following sequence to test the 
convergence of the basis set. Starting with the simplest 
scheme, a single Q (SZ) basis, a second group of valence 
orbitals is added for flexibility, forming the double-^ (DZ) 
basis. Por completeness, we also add polarization orbitals 
to both valence sets, generating SZP and DZP bases. 

Finally, it is possible to optimize the localization radius 
in order to increase accuracy for a given basis set. While 
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TABLE I: Comparison between converged basic parameters 
for bulk GaAs. a, B and Etotai represent the lattice constant 
(A), the bulk modulus (GPa), and the total energy per GaAs 
pair (eV/pair), respectively. Eg denotes the energy gap (eV) 
and A_ff is the heat of formation of GaAs calculated using 
Eq. 2](eV). See the text for description of the basis set la- 
bel. Results are compared to a recent plane-wave calculation 
(PW)^^ and to experimental values from Ref'^ at OK, unless 
other references are cited. 





SZ 


SZ-0 


SZP-O 


DZ 


DZP 


PW 


Expt. 


a 


5.68 


5.66 


5.6 


5.64 


5.6 


5.55 


5.65 


B 


59.3 


68.9 


78.8 


67.7 


70.4 




75.3 


Etotai 


235.2 


235.5 


235.9 


235.7 


236.0 








0.61 


0.78 


0.98 


0.66 


0.82 


1.08" 


1.54 














0.7'' 




AH 


0.66 


0.99 


0.78 


0.81 


0.72 


1.0" 


0.73= 



0.67^^ 
0.83^* 

" Reference^ '' Reference Reference 
Reference 12 " Reference 13 



the computational efficiency is slightly reduced, as the 
optimized orbitals have generally a longer tail, it is often 
a good alternative t o in creasing the size of the basis set 
(for details see Ref. 1^30 ). We also test the accuracy of 
these basis sets optimized with respect to the amount of 
overlap between atomic orbitals around the defect using 
the optimizing procedure of Anglada et al.^ at 0.0 GPa 
for the SZ basis set. We find that the efficiency of these 
orbitals with (SZP-O) and with (SZ-0) is comparable to 
those of DZP and DZ, respectively. 

Table Q reports the values of a number of structural 
and thermodynamical quantities for bulk GaAs as com- 
puted using these various bases with a /c-point sampling 
density of O.OSA"^, corresponding to that for a 216-atom 
unit cell with a 2 x 2 x 2 fc-point sampling. For SZ, the 
lattice constant at zero pressure is found to be 5.68 A, 
overestimating the experimental value by only 0.03 A. 
The density increases with the size of the basis set and 
the lattice constant for DZP is found to be too small by 
0.05 A with respect to experiment. The relatively con- 
tracted structure obtained with DZP is characteristic of 
LDA; plane-wave calculations also using LDA give 5.55 
Increasing the basis set leads to a significant improve- 
ment on the calculated value of the bulk modulus as it 
goes from 59.85 GPa for SZ to 70.4 GPa for DZP, close 
to the experimental value of 75.3 GPa. 

The LDA band gap is found to be 0.61 eV, 0.66 eV and 
0.82 eV for SZ, DZ and DZP basis sets respectively, un- 
derestimating, as usual with with this approximation, the 
experimental gap of 1.54 eV. The DZP band gap lies well 
within the range of energy gaps obtained from PW (0.7- 
1.0 eV)f2ii£ however, and can be considered converged. 

No such systematic problem is found for the total en- 



ergy and the heat of formation. In particular, the heat of 
formation obtained with DZ and DZP is very close to the 
experimental valuei^, showing a better agreement than 
previous plane-waves calculations&iiii^ 

Overall, therefore, we see a well-defined trend in the 
structural and thermodynamical values shown in TableQ] 
most quantities converge rapidly as the basis goes from 
SZ to DZ to DZP, with DZP providing an excellent agree- 
ment with experiment. Moreover, it appears that the 
optimal basis sets, SZ-O and SZP-O, compare very well 
with DZP, suggesting that they could be used when com- 
putational costs are an issue. The application of these op- 
timized bases to study the diffusion of vacancies in GaAs 
will be reported somewhere else^^ 

B. Defects formation energies in supercell 
calculations 

When computing structural and energetic properties of 
defects using ab-initio methods, it is important to ensure 
that the size of the basis set is complete enough but also 
that the simulation cell is sufficiently large to avoid self- 
interaction between the defect and its images. We have 
shown previously, in a study of the neutral vacancy in 
silicon, that a supercell of at least 216 atomic sites can be 
necessary in order to reduce the elastic and electronic self- 
interaction and obtain the right symmetry around the 
defect»iSi As discussed in Sect. IIVB II we find a similar 
behavior for As vacancy; unless indicated, therefore, we 
use a 215-atom cell for all our calculations of defects. 

In all calculations, this initial 215-atom configuration 
cell is randomly distorted, to avoid imposing spurious 
symmetry in the fully relaxed defect state. All atoms are 
allowed to relax without any constraint until every force 
component falls below 0.04 eV/A. The energy minimiza- 
tion takes place at a constant volume, using the optimal 
lattice constant obtained with DZP, 5.6A (see Tab. QJ, 
1% denser than the experimental value. 

The formation energy can be evaluated directly from 
total energies obtained from electronic structure calcu- 
lations. For binary compounds it is current to use the 
formalism of Zhang and Northrupi2i (see Ref. for in- 
termediate steps). The formation energy of a defect of 
charge state q is defined as: 

^ + ^(^v + Me) - ^{nAs - ^100)^1^- (1) 
where E j is independent of A/i and /.te, and is represented 

by 

-linA.-nGa){^^T~^^Gl') (2) 

where uas and are the number of As and Ga ions 
present in the sample, q denotes the net number of elec- 
trons or holes supported by the vacancy, Me is the electron 



3 



chemical potential or the Fermi energy Ep, and Ey is 
the energy at the valence band maximum. Errors in Ey 
due to the finite supercell are corrected by aligning the 
vacuum levels of the defective supercell and undefected 
superccll^^. 

If A/i is defined as the chemical potential difference: 

Am = {^iAs ~ ma) ~ inT - y^Ga ). (3) 

the restriction on the chemical potentials becomes: < 
Me < Eg and —AH < Am < AH, where Eg is the en- 
ergy gap, and the heat of formation AH of bulk GaAs is 
defined as the difference between the chemical potential 
of bulk As and bulk Ga crystals and that of bulk GaAs. 
This latter quantity represents the energy necessary to 
dissociate GaAs into its individual components: 

A C7 ..bulk , bulk bulk /a\ 

AH - Mas + MCa - MCaAs (4) 

For Ga vacancies of charge q the Eq.^ reduces to : 

Ef = e) + q{Ev + Me) + ^Am (5) 

with E) ^ EUq) - '^P^GaAs + - y^GaY for As 

vacancies, it becomes: 

Ef^E!^+q(Ey+p.,)-\^A^l (6) 

and E^ = E^^M) - ¥/^Gal " Wa^ " M^^f')- 



±3, respectively. The quadrupole correction, which we 
evaluated by hand, is proportional to Xjl? . For the 215- 
atom supercell, it is 2.42 x 10^^ * qQ eV (where Q is the 
quadrupole moment), and can therefore be neglected (see 
also Ref.m. 

Because of the limitations of LDA, localized DFT 
eigenvalues are not equivalent to the measured electronic 
levels. Thus ionization energy is obtained from the differ- 
ence between q\ and 02 electron total energy calculations 
(e(g2/gi) = EH^-Ef^^-{q2-q\)Ev), rather than the dif- 
ference of qi and q\ electron eigenvalues of a single calcu- 
lation. Usually only one electron is transferred between 
the electron reservoir and the defect levels. When two 
electrons are transferred at the same time the electron- 
electron repulsion is compensated by a relaxation of the 
structure around the defect that arises from a strong 
electron-phonon coupling. This so-called negative-U ef- 
fect is found when the ionization level e{q — f/q) appears 
above £{q/q+l), thus a direct transition (q—l/q+l) is 
energetically more favorableiS^ 



III. CONVERGING DEFECT FORMATION 
AND IONIZATION ENERGIES 

In this section, we study the effects of the basis set and 
the /c-sampling and the simulation cell on the formation 
energy and the relaxed geometry of neutral and charged 
Ga and As vacancies. 



C. Computing the ionization energy of charged 
defects 

The concentration of charged defects is controlled by 
the position of the Fermi level which is determined by 
the local concentration of carriers. Since GaAs is used 
in a doped state in devices, it is important to assess the 
possible charged states of defects. 

Charges can affect strongly the formation energy as 
well as the structure of a defect, changing the symme- 
try of the relaxed state and altering considerably the lo- 
cal electronic properties. For charged defects, the effects 
of finite-size supercell will be even more marked due to 
long-ranged nature of the Coulomb interaction; the use a 
sufficiently large supercell is therefore even more impor- 
tant. 

To account for the electrostatic interaction of period- 
ically arranged defects of charge q as well as their in- 
teraction with the compensating background, we follow 
the approximate procedure of Makov and Payne. The 
correction to the total energy of a charged system is han- 
dled by SIESTA, and it consists of a monopole correction 
only ( q^a/2eL), where a is the Madelung constant of 
the simple cubic lattice, L is the defect-defect distance 
(16.8 A) and e is the experimental static dielectric con- 
stant. The monopole correction is found to be 0.094 eV, 
0.37 eV and 0.84 eV for the charge states ±1, ±2, and 



A. Local basis set effect 

We first study the effect of the choice of the local basis 
set in the defects. In order to cancel out all other effects, 
we use supercells of 215 atomic sites and a Mokhorst- 
Pack grid^"^ of 2x2x2 corresponding to a density of 
0.03A~i of k-points. 

Cohesive energies as well as bulk moduli studied in 
Section III Al for the different basis sets used show that 
atomic bonding is strengthened progressively as we go 
from SZ basis to DZP. 

Structural relaxation is directly related to the inter- 
atomic forces acting on the atoms around the vacancy 
and on the strength of atomic bonding. Atoms around 
the vacancy form initially an ideal tetrahedron with six 
equal distances labeled di — with tetrahedral symme- 
try Td- After the full relaxation (see Sect. Ill All , distances 
and angles can be altered and the symmetry is either con- 
served or broken. 

Relaxations around the vacancies are given for different 
charge states in Tables Hll and Hill for Ga and As vacan- 
cies respectively. In both cases the corresponding forma- 
tion energy is reported as well as the relative change in 
the volume of the tetrahedron with respect to the ideal 
oneii^^ 

Due to the finite precision in the relaxation, there is 
some imprecision in the identification of the defect sym- 
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metry. Here, if the highest relative difference between 
two bonds is lower than 1 % (equivalent to a precision 
of 0.04 A), the structure is assigned to highest symmetry 
group. The last column lists the symmetry groups for the 
different defects in the DZP basis set. Unless specified, 
the symmetry group for all bases is the same as that of 
DZP. 



1. Vaa 

The Ga vacancy maintains the same symmetry for all 
charge states irrespective of the basis set used. The struc- 
tural relaxation is most important for the smallest basis 
set, SZ, decreasing progressively by about one third as 
the number of bases is increased, but the Td symmetry is 
maintained in all cases, with the atoms moving inwards 
systematically. Moreover, for each basis set the degree of 
structural relaxation, is almost independent of the charge 
state, contracting very slightly, by about 2 %, from the 
neutral to the -3 charge state, for all basis sets. 

As can be seen in Tabled the formation energy is also 
rather well converged with the minimal basis set (SZ), by 
comparison with the more accurate DZP: the difference 
between the two bases is at most 0.2 eV. 

Heavily charged defects like Vq^^ which is the most 
likely charge state in a heavily doped material, are more 
sensitive to the completeness of the basis set. This ef- 
fect explains the fact that the formation energy decreases 
with an improved basis, contrary to the other defects. 

Because of its technological importance, we must en- 
sure that the orbital overlap around the vacant site is 
sufficient to accommodate the extra electrons in Vq^. 
We can do so by placing a ghost Ga atom at the defect 
site. For that purpose, we generate a set of orbitals and 
place them on the crystalline site, without adding any 
pseudo-potential or extra electrons. The system is then 
relaxed using the same convergence criterion as before. 
This ghost atom does not have any significant effect on 
the total energy of the defect when the DZ and the DZP 
basis sets are used. In contrast, SZ basis total energies 
are corrected by 0.34 eV, suggesting that the SZ orbital 
are too short. This additional set of orbitals is sufficient 
to correct for the overestimated formation energy with 
SZ basis, decreasing its value from 4.06 to 3.72 eV, fol- 
lowing the general trend observed for other charge states 
(see table HH). 

2. Vas 

The situation is very different for the As vacancy: the 
local symmetry is broken for most charge states and the 
completeness of the basis set impacts strongly on the re- 
construction around the defect. Except for the positively 
charged vacancy, the bonds are stretched considerably to 
form pairs, leading to volume deformation by as much as 
60%. Because of this strong deformation, we must relax 



the criterion on the symmetry. To allow the reader to 
judge the classification, we write the distances di — in 
ascendent order in the last column of table IIIII 

As a general trend, the SZ basis is less efficient for As 
than for Ga vacancies: the formation energy is underes- 
timated by as much as 0.6 eV (27%) for charge +1 and 
1.0 eV (17%) for charges -3 as compared with DZP. The 
underestimation drops to 2-8 % with DZ, a considerable 
improvement, for all charge states. The improvement 
in formation energy can be directly correlated with the 
prediction quality of the local relaxation. For example, 
while Td symmetry is conserved for the three basis sets for 
V^^ , the change in the volume around the V^/ is highly 
overestimated by 135 % using SZ basis compared to the 
DZP results. The overestimation drop to 35 % with DZ, 
leading to an error of less than 0.1 eV compared with 
DZP. 

As more electrons are added to the defect level this 
trend tends to diminish, the symmetry and the relaxation 
of the defect can be described with reasonable accuracy 
using the DZ basis. Applied to the singly negative As 
vacancy V^/, the inclusion of the ghost atom at the va- 
cant site has a smaller impact on the energy level and 
relaxation using SZ than for Vq^ , and next to none with 
DZ and DZP. For the minimal basis set, the correction is 
negligible and accounts only for 0.067 eV. 



3. Convergence of the ionization energy with the basis sets 

As defined in section lTl Gl ionization energies are taken 
as the difference between the total energy in different 
charge states. Unlike formation energies: the errors in- 
troduced for the ionization energies are considerably re- 
duced errors coming from LDA are canceled out and 
those coming from different chemical potentials are elim- 
inated. The only remaining error are those coming from 
the basis set convergence, the evaluation of the valence 
band edge energy and the Madelung correction. 

Figure n displays the convergence of the ionization en- 
ergies as a function of the basis set used for Vas and 
Vca- For both types of defects, the preliminary results 
obtained using SZ basis set give a rough estimate of the 
location of ionization energies in the band gap. These 
energies converge with increasing basis sets, but slower 
than the formation energies. For Vca, while the forma- 
tion energies are reasonably converged already with the 
minimal basis set, the ionization levels found using the 
SZ basis are noticeably overestimated compared to DZP. 
The levels are located at 0.21, 0.43 and 0.69 eV for SZ, 
while for DZ it changes for 0.13, 0.35 and 0.57 eV then 
finally converges to 0.05, 0.4 and 0.55 eV for DZP. For 
Vas, we found in section IlII A 21 that formation energies 
are underestimated by less than leV with SZ for most 
charge states. Once ionization energies are calculated, 
errors on the formation energies coming from the relax- 
ation around the defect cancel out since defects in all 
charge states suffer from this effect. 
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TABLE 11: Convergence of the formation energy in eV with respect to the basis set for the Ga vacancy. Relaxation around 
the vacancy are given in % compared to the ideal tetrahedral distance between As nearest neighbors. The distances are labeled 
di—de, the negative sign indicates an inward relaxation. The tetrahedron volume change is also given in % of the ideal volume 
{AV = 100 * {V — Vo)/Vo)- The last column displays the symmetry group of the defect (see the text for more details). 



Distances in % 



Basis 




di 




d-i 


di 


ds 


de 


AV 


Symmetry 






















SZ 


2.7 


-19.3 


-19.2 


-19.3 


-19.3 


-19.2 


-19.3 


-47.3 




DZ 


2.8 


-14.5 


-UA 


-14.5 


-14.4 


-14.5 


-UA 


-37.4 




DZP 


2.9 


-13.5 


-13.4 


-13.5 


-13.4 


-13.4 


-13.5 


-35.2 


Td 












T 1 










SZ 


2.9 


-19.0 


-19.1 


-19.0 


-19.0 


-19.0 


-19.0 


-46.9 




DZ 


2.9 


-14.7 


-14.7 


-14.7 


-14.7 


-14.7 


-14.7 


-37.9 




DZP 


3.0 


-14.2 


-14.2 


-14.2 


-14.2 


-14.2 


-14.2 


-36.9 


Td 






















SZ 


3.4 


-19.5 


-19.5 


-19.5 


-19.5 


-19.5 


-19.5 


-47.8 




DZ 


3.2 


-14.4 


-UA 


-UA 


-UA 


-14.4 


-14.4 


-37.3 




DZP 


3.4 


-14.0 


-14.1 


-14.0 


-14.0 


-14.1 


-14.0 


-36.5 


Td 












''Ga 










SZ 


4.1 


-19.7 


-19.8 


-19.7 


-19.7 


-19.8 


-19.7 


-48.3 




DZ 


3.8 


-15.1 


-15.2 


-15.1 


-15.1 


-15.2 


-15.1 


-38.9 




DZP 


3.9 


-14.6 


-14.6 


-14.6 


-14.6 


-14.7 


-14.6 


-37.7 


Td 



SZ basis gives a correct qualitative description of the 
nature of the electronic transition (double negative-U ef- 
fect, to be discussed in section llVB The location of 
the levels in the band gap as well as the distance be- 
tween them are also reasonably converged. Ionization 
energies with SZ are slightly underestimated compared 
to the DZP ionization levels. Levels are at 0.13 and 1.19 
eV for SZ , 0.24 and 1.09 for DZ and finally 0.27 and 1.27 
for DZP. 



By studying the effect of the choice of the basis set 
on the structural relaxation, the formation and the ion- 
ization energies, we conclude that the SZ basis is signif- 
icantly less efficient when the local symmetry is broken. 
In these cases, the use if a second radial function (DZ) 
is necessary to obtain reasonable numbers. Moreover, 
both DZ and DZP are complete enough to represent the 
properties associated with a defect without the need for 
a ghost atom. 



B. K-points effects 

The effects of Brillouin zone sampling are studied by 
comparing the formation energies and the relaxation vol- 
ume of both As and Ga vacancies for all charge states. 
We consider two k-point sampling, a F point sampling 
and a 2 X 2 X 2 Monkhorst-Pack mcsh^"^ , corresponding, 
for the 215-atom cell used here, to a density of O.Q&K~^ 
and O.OSA"^, respectively. To isolate sampling effects, all 
calculations are done using the DZP basis set. Results 
are shown in Table HVl 

The use of T point only gives formation energies that 
are reasonably converged for Ga vacancies. Similarly, 
the bond lengths around the vacancy and the relaxation 
volume (Fig. 121) are near those of the 2x2x2, preserving 
the Td relaxation symmetry for all charge states. 

The use of F-point only produces less reliable results 
in the case of As vacancies. Formations energies are well 
converged but the relaxation symmetry around the defect 
is not correctly predicted for all charge states. With F- 
point, the distance between the atoms forming the tetra- 
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TABLE III: Convergence of the formation energy Ej, in eV with respect to the basis set for the As vacancy. The relaxation 
around the vacancy is given in % compared to the ideal tetrahedral distance between Ga nearest neighbors. The distances are 
labeled bi — &6, the negative sign indicates an inward relaxation. The volume change around the vacancy is also given in % of 
the ideal volume {AV — 100 * (V — Vo)/Vo). The last column displays the symmetry group of the defect (see the text for more 
details) . 



Distances in % 



Basis 


Jif[e\ ) 




0,2 


"3 


0,4 


"5 






Symmetry 






















sz 


2.2 


-14.5 


-14.9 


-15.1 


-15.7 


-15.9 


-16.6 


-39.5 




JJZj 


Z. I 


- / .y 


-o.U 


Q 1 
-O.i 


-0.4 


-0.4 


-o. / 


-zz. / 




J_/ Zj A 




-O.O 


-o.y 


fl 1 
-0. 1 


-o.o 


ft o 


ft A 


17 1 


-Id 






















sz 


2.4 


-13.8 


-14.3 


-15.2 


-15.2 


-31.3 


-31.7 


-53.4 


Cih 


T\'7 


O 1 
O.i 






1 O '7 


1 O A 

-io.U 


-ZD.y 


-z 


-4d. / 








11 1 


1 1 O 


- J. J..O 


1 1 R 
- J. J..O 


0/1 ft 




























SZ 


2.5 


-14.3 


-14.7 


-15.3 


-15.4 


-33.3 


-33.5 


-55.6 




T\'7 


6.2 


IOC 




1 o o 

-io.z 




-oi.o 


-oz.U 


KO /I 

-Oz.4 




LJZuir 


o.o 


lift 


1 1 Q 




1 O ft 


Q1 *> 


Q1 PI 


Pll o 


~ ^2d 






















sz 


3.8 


-17.8 


-18.5 


-18.5 


-19.3 


-33.6 


-33.7 


-58.6 




DZ 


4.4 


-16.5 


-17.1 


-17.5 


-18.2 


-32.0 


-32.3 


-56.2 


Csi, 


DZP 


4.6 


-16.5 


-16.8 


-17.0 


-17.6 


-31.6 


-31.7 


-55.4 


C:iv/D2d 






















SZ 


4.9 


-18.4 


-18.4 


-31.4 


-31.6 


-31.7 


-31.8 


-65.5 




DZ 


5.4 


-18.3 


-18.2 


-30.4 


-30.6 


-30.7 


-31.0 


-63.8 




DZP 


5.9 


-18.1 


-18.2 


-30.4 


-30.2 


-30.2 


-30.2 


-62.9 


i32d-resonant 



hedron as well as its associated volume are already con- 
verged for charge stated -1-1, 0, and -1, as shown in Fig.|21 

For , F-point sampling overestimates the volume 
contraction around the defect by 7.5%. The long bonds 
in the pairing mode associated with the D2d structure are 
shortened from 3.46 A for this sampling to 3.26 A for 2 x 
2x2 sampling. The k-point sampling effect becomes even 
stronger for the highly charged V^^. The resonant bond 
geometry present with high-density sampling (discussed 
further in Sect. IIVB ip is not found with the F-point 
sampling and the defects relaxes into the usual pairing 
mode with two short bonds and four long bonds. 

Comparing the two samplings for V^^, we find that: 
(1) with DZP + 2^ mesh, the resonant-bond configura- 



tion has the lowest energy; (2) the pairing mode con- 
figuration is unstable using DZP + 2^ mesh, starting 
from this configuration, the defect relaxes back into the 
resonant-bond state; (3) for F-point sampling the reso- 
nant bond configuration is found to be metastable but 
the pairing mode is favored with an energy difference of 
0.6 eV. 



Overall, the use of F-point sampling gives a reasonable 
description of the energetic and structural properties of 
the Ga vacancy as well as most charge states for As. It 
fails, however, for heavily charged defects, when many 
electrons are involved in the bonding. 
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FIG. 1: Schematic representation of the convergence of ion- 
ization energies as function of the basis set. Due to the un- 
derestimation of the gap, as a consequence of LDA, it is usual 
to align the conduction band maximum (CBM) with the ex- 
perimental value, at 1.54 eV from the valence band maxi- 
mum (VBM). The left panel displays the ionization energies 
as function of the basis set for the Ga vacancy that are all 
all below the experimental mid gap. The right panel shows 
the two ionization levels located below and above the mid gap 
and their convergence with respect to the different basis sets 
used. Refer to the text for the basis set description. 



TABLE IV: Convergence of the formation energies (in eV) 
for Ga and As vacancies with respect to the Brillouin zone 
sampling. 





+ 1 


Neutral 


-1 


-2 


-3 








EfiVoa) 
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2.51 


2.73 


3.17 


3.81 


2x2x2 




2.94 


3.00 


3.40 


3.94 






15% 


9% 


7% 


3% 














r 


2.73 


3.21 


3.52 


4.70 


6.08 


2x2x2 


2.79 


3.25 


3.33 


4.52 


5.86 




-2% 


-1% 


6% 


4% 


4% 



C. Size effects 

Size effects on defect formation energies liave been 
widely discussed for a number of systems including the 
silicon vacancy (see Ref. and references therein) and 
GaAsi^ Size effects are found to be strong for cubic su- 
percells smaller than 216±1 for both materials and be- 
comes negligible for larger systems. For this reason, we 
will study in this section size effects on charged Ga and 
As vacancies for supercells of 216 ionic sites and smaller. 



FIG. 2: (Color online) Histogram of the change in the volume 
of the relaxed tetrahedron formed by atoms surrounding the 
vacancy in % of the ideal one for two densities of k-points in 
the Brillouin zone (refer to the text for details). For compar- 
ison, data for Ga and As vacancies in different charge states 
are plotted in the same figure. 



To characterize the size effects, we consider the domi- 
nant charge state for each vacancy type, V^^ and Vq^, 
with DZP basis set. We simulate two cell-size with the 
same k-point density: a 63-atom unit cell with a 3x3x3 
sampling and a 215-atom cell with 2x2x2 sampling iS4 

We find a formation energy of 3.83 eV compared to 
3.33 eV obtained for the 215 cell for V^/. In this case 
size effects are strong and overestimate the formation en- 
ergy by 0.5 eV. More important, the symmetry of the 
relaxed defect is different for the two sizes: in the 63- 
atom cell, the defects relaxes to a C21, symmetry while 
the defects 215-atom cell adopts a Z?2d symmetry (Ta- 
ble IIII|I . We check that the 63-atom cell is not caught 
into a metastable state by relaxing the defect starting in 
a D2d symmetric state. After full relaxation, the cell re- 
laxes into a C2V , confirming that this is the lowest-energy 
symmetry for this cell. The different symmetry also im- 
pacts the change in the volume surrounding the defect. 
As a consequence of the strong defect-defect interaction 
in the smaller cell, the structural relaxation is hindered 
and the volume decreases by 43.8% compared with 51.2% 
for the 215-atom cell. 

As could be expected from the previous sections, size 
effects are less important for Vq^. In particular, the ini- 
tial tetrahedral symmetry (T^), with atoms equidistant, 
is maintained around the defect for the relaxed 215-atom 
supercell. The same symmetry is found in the 63-atom 
supercell. Intercell defect-defect interaction rigidity the 
lattice, however, and the change is volume is only 34.8 
% for the 63-atom cell compared to 37.7 % for the larger 
supercell. As with V^/, the formation energy is over- 
estimated with the small cell: Ej^ — 4.24 eV while 
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= 3.94 eV). This difference is considerable as it 
is on the order of the ionization energies. 

In this section, we have studied in details the effects of 
the choice of basis set, k-point sampling and simulation 
cell on the properties of charged defects. In summary, we 
find that: (1) The DZP basis set is well-converged and 
ensures reliable results for all charge states. For a num- 
ber of charged states, it is also possible to use cheaper 
optimized basis sets for a similar accuracy. This is not 
always the case, however, and the applicability of these 
basis sets must be evaluated on a case by case basis. (2) 
SZ is less efficient than DZP for defects where symmetry 
is broken, but it gives a satisfactory estimation of the 
location of ionization levels in the band gap. (3) For su- 
percell of 216 atoms or more, the density of k-points has 
only a minor effect on the defect relaxation. The use of 
the r point only gives a relaxation and a symmetry that 
are satisfactory, in most cases. However, this reduced 
sampling must be used with care for highly charged de- 
fects such as V^^ . (4) The errors arising form size ef- 
fects are much more important than those coming from 
the density of k-points. In particular, size effects can be 
the source of errors in estimating the ionization energies, 
especially when the transition from a charge state to an- 
other induces breaking of the symmetry. 



TABLE V: Comparison between ionization energies (mea- 
sured from the valence band edge) of the Ga vacancy in GaAs. 
Results are grouped following the three classes discussed in 
the text. 



Authors 




Ionization levels (eV) 


U/-1 


-1/-2 


-2/-3 


Negative U 
+ 1/-1 -1/-3 


Seong and Lewis— ^ 








035 078 


Northrup and Zhang2£ 


0.19 


0.2 


0.32 




Poykko et a& 


0.11 


0.22 


0.33 




Schick et ali^ 


0.09 


0.13 


0.2 




Janoti et al. — - 


0.13 


0.15 


0.18 




Jansen and Sankey— - 


0.1 


0.35 


0.50 




Baraff and Schliitej— 


0.2 


0.5 


0.7 




Cheong and Chang— 




0.49 


0.69 




Gorczyca et al— 


0.39 


0.52 


0.78 




This work 


0.05 


0.4 


0.55 





IV. RESULTS AND DISCUSSION 

We discuss here the results reported in the previous 
sections; we concentrate on the highly converged results 
for the Ga and As vacancies, obtained using the DZP ba- ^ 
sis set using a 215-supercell and a Monkhorst-pack grid 
of 2 X 2 X 2 in the reciprocal space. We invite the reader ^ 
to refer to the DZP results reported, in bold, in Tables UTI u 
and mil We first deal with the stochiometric case where S 
nca = nAs (A/Lt =0), then we study the formation ener- S 
gies for each kind of defects separately under ideal grow- '-^ 
ing conditions as a function of the doping level. Next, g 
the dominant vacancy type defects in real GaAs crystals o 
are identified by taking into account growing conditions. 
The ionic chemical potentials varies from As-rich condi- 
tions (A/i — —AH) to Ga-rich conditions (A/^ = +AH) 
as the Fermi level is changed progressively. 



A. Gallium vacancies 
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1. Relaxation geometry 



Table describes the fully relaxed geometry of the 
defect. The structural deformation obtained is well lo- 
calized around the vacant site; the magnitude of the re- 
laxations is listed for the nearest neighbor As atoms in 
all relevant charge states, reaching 15% of the bulk bond 
distances. The tetrahedral symmetry Td is always con- 
served for this defect, irrespective of the charge state: 
only the breathing mode matters here. The As dangling 



FIG. 3: (Color online) Formation energies as function of Fermi 
level in various charge states of Ga vacancies at OK. The Fermi 
level is defined by reference to the valence band maximum. 
Ionization levels are defined as the intersection between the 
formation energies of different defects. Defect with the lowest 
formation energy is dominant. Arrows point to the location 
of the ionization levels labeled (a) for (0/-), (b) for (-/-2) and 
(c) for (-2/-3) 
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bonds do not form pairs in any charge state, but the 
back bonds formed with Ga atoms are clearly weakened; 
this is in agreement with the observation that the pairing 
mode is generically not energetically favorable for cation 
vacancies (Ga),"^^ We also observe that all As atoms re- 
lax inward but the amount of the relaxation does not 
increase as more electrons are added to the vacancy lev- 
els but remains stable (13.5-14.5%). Our results are in 
good agreement with the results of Laasonen et aliiS, 
and Seong and Lewia^, who found a systematic inward 
relaxation with tetrahedral symmetry for Ga vacancies 
in (0,-1,-2) charge state studied using an empirical tight- 
binding potential. 



2. Energetics 

Experimentally, the Ga vacancy is found to exist in the 
0, — ,— 2,— 3 charge states. Until recently, the preferred 
charge state for the Ga vacancy in GaAs was the subject 
of a hot debate: most ab-initio calculationsiSiSiSSiSl find 
that where GaAs in either semi-insulating or n-type and 
Fermi energy is away from the valence band edge, the 
gallium vacancy is in the triply negative charge, while 
diffusion experiments suggest a charge of —2 or —1.^^ 

Using positron annihilation to determine the Gibbs 
free energy of formation for Ga vacancies in GaAs, 
Gebauer et alm^ could finally resolve this debate, giv- 
ing a quantitative estimation of the formation enthalpy 
for Vq^ and Vq^. The vacancy concentration is directly 
probed with positron annihilation in Te doped GaAs as 
function of doping concentration, temperature and chem- 
ical potential. Our estimate of the formation energies re- 
ported in the first column of Table ^ are in good agree- 
ment with recent experimental and theoretical data. For 
the neutral vacancy, we get Ej{Vq^)= 2.94 eV, a val- 
ues that agrees with E ^{Vq^)—2.8 eV from Bockstedte 
and SchefHer— (earlier first-principals calculations^^ pre- 
dict 3.5 eV) and the experimental results of Gebauer et 
al.^ H'f{V^J=3.2± 0.5 eV. A formation enthalpy of 
1.8±0.5 eV was also measured by Mitev et alm^ using 
inter-diffusion experiment on AlGaAs/GaAs heterostruc- 
ture. However, the charge state of the associated defect 
is unknown, rendering the comparison with our results 
difhcult. 

For the triply negative charge state, we find that 
e'j:{V^^) = 3.9 eV. This value compares well with re- 
cent experiments-'^'^ which give E[j{Vq^) = 3.6 eV as well 
as with recent theoretical study by Janoti et al. — ( 3.6 
± 0.2 eV ). 

Considering the stochiometric case where nca = nAs 
(A/x =0), we can study the formation energies of GaAs 
under ideal growing conditions as a function of the doping 
level. FigureOdisplays the formation energies as function 
of Fermi energy for various charge states of Ga vacancies 
at OK. As Ef depends linearly on the electronic chemical 
potential, /ie, the slope of each of the lines represents the 



TABLE VI: Comparison between ionization energies (mea- 
sured from the valence band edge) for the As vacancy in 
GaAs. The data are grouped according to the two categories 
discussed in the text. 



Authors 


Ionization levels (eV) 


+2/+1 +1/0 


Negative U 
0/-1 +1/-1 -1/-3 


Seong and Lewis— ^ 


1.41 


1.54 


Jansen and Sankey— - 


1.3 


1.40 


Cheong and Changi^ 




0.785 


Poykko et al2i 




0.86 


This work 




0.27 1.27 



net charge for the system q. Intersections determine the 
ionization levels where one electron is transferred from 
the electron reservoir to the defect level. We can see that 
at each transition only one electron is transferred at a 
time. Moreover, the ionization levels labeled a, b, and 
c favor the stability of the —3 charge states for intrinsic 
and n-type GaAs. 

In order to compare with earlier theoretical work, 
we summarize the results into three types of behavior: 
(1) Sole among all calculations, Seong et Lewis^^ find 
a negative-U effect for the Ga vacancy using TB-MD 
method, the transition levels identified are very shallow 
and favor the triply negative state in almost the entire 
range of the Fermi l evel. (2) Th e second category of lev- 
els reported in Ref. Illl24l26l27l are shallow and lie well 
below the mid gap with no negative-U effect detected. 
(3) Defects level can also lie deeper below the mid gap. 
This is the case for levels a, b and c in Fig. O which 
compare well with the results of Cheong and Chang^^, 
Baraff and Schliiter— and a recent study form Gorczyca 
et al"^*^ where the — 3w charge of Vca is relevant only if 
the Fermi energy is above 0.55 eV. 

Unfortunately, ionization levels cannot be directly 
measured experimentally, only their sum is obtainable. 
Nevertheless, experimental values obtained by electron 
irradiation of GaAs^^ support the assignment of deep- 
lying levels in GaAs, in agreement with the third cate- 
gory. Gebauer et al.^^ confirm this assignment using a 
model to fit their experimental data in order to identify 
the charge state of the vacancy in GaAs from the loca- 
tion of the ionization levels. This calculation shows that 
the -3 charge state is the most stable charge state if deep 
ionization energies coming from a calculation on unre- 
laxed Ga vacancies of Baraff and Schliiterj^S which are 
in agreement with our more precise calculations. 
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FIG. 4: (Color online) formation energies as a function of 
the Fermi energy of various charge states of As vacancies at 
OK. The Fermi level is calculated with respect to the valence 
band maximum. Arrows point to the location of the ionization 
levels labeled (a) for (+/-), (b) for (-/ - 3) 



B. Arsenic vacancies 



1. Relaxation geometry 

As vacancies in various charge states have been stud- 
ied by a number of authorsi2i2Si2i taking into considera- 
tion ionic relaxation and relaxation geometry. Except for 
^As' breathing mode displacement breaks the vacancy 
local symmetryiifliS^ Contrary to Ga vacancy, the vol- 
ume of the As vacancy increases as electrons are added 
and shrinks when electrons are removed. 

Our results confirm part of these findings. From Ta- 
ble mil the volume of the tetrahedron shrinks from -17% 
to -60% of the initial volume as electron are added to the 
vacancy. However, all charge states, even positive ones, 
display inward relaxation with respect to the unrelaxed 
volume. 

Comparing the bond lengths obtained after a full re- 
laxation of the structure (Table [ffl}, we find that for the 
positively charged vacancy (V^/) there is no electron in 
the localized states and all four atoms relax inward by 
about -6% conserving the tetrahedral symmetry, with no 
Jahn- Teller distortion. Although this inwards relaxation 
is more important, in absolute value, than the outward 
relaxation reported in referencesAfl^^^, they agree in term 



of the conservation of the symmetry. 

For the neutral As vacancy in GaAs {V^s), there is just 
one electron in the localized state formed by the dangling 
bonds. The volume reduction is more than twice as large 
as for the positively charged defect. This change in vol- 
ume is associated with a Jahn- Teller distortion with the 
formation of two dimers, leading to two short (-11.5%) 
and four long Ga-Ga bonds (-24.5%) arranged in a D2d 
symmetry. This stretches and weakens the back bonds 
but it allows all atoms to recover a fourfold coordination. 

Using ab initio molecular dynamics, Laasonen et al.^ 
see a small (2-3%) outward relaxation, and an even 
smaller (0.6%) pairing-mode relaxation, leading to a 
weak tetragonal distortion with D2d symmetry. This cal- 
culation was found to suffer from band dispersion for the 
localized defect states due to the artificial interaction 
between unit cells; as a consequence atoms surround- 
ing the vacancy are not allowed to relax properly. Us- 
ing tight-binding molecular dynamics, a larger breathing- 
mode displacement was obtained by Seong and Lewis^ 
They find that local tetrahedral symmetry was broken, 
as one neighbor atom of the defect relaxes inward while 
the other three relax outward. The pairing mode relax- 
ation is also found to be very small. Feng et al.^ using 
a similar method with a 64-atom supercell, found similar 
results leading to trigonal distortion with C^y symmetry. 

The difference between our results and previous calcu- 
lations come from the use of a better converged poten- 
tial as well as of a larger unit cell. Moreover, as stated 
previously, we have started the relaxation from various 
random geometries, always converging to the same final 
state: the lowest energy configuration has D2d symmetry 
for neutral Vas- 

The extra electron added to get V^^ can be accom- 
modated in the same localized level as the previous one. 
A stronger pairing mode relaxation appears and the two 
short bonds become stronger (from —24.5% to -31.5%) 
atoms of the dimer get closer, while the long bond are al- 
most kept fixed. In agreement with our results reported 
in Table IITTI Chadi,^''- using LDA and 32-atom supercell, 
finds that the -1 charge state arises from a direct transi- 
tion from the +1 state due a pairing of the neighboring 
Ga atoms and has C2v symmetry (pairing mode). The 
Ga atoms would then move by 0.8 A (- -20%) from their 
ideal position to form two sets of paired bonds. 

The next electron added to the localized levels, occupy 
a different state. The arrangement of the atoms around 
the doubly negative vacancy Vj^^ is directly affected: the 
short bonds remain unchanged but the long bond be- 
come stronger passing from —12.6% to —17.6% with a 
slight change in the relaxation volume. The dimers are 
therefore brought closer without affecting the intradimer 
distance. 

Most interestingly, the relaxation geometry is modified 
when a fourth electron is added. In particular, there is 
an inversion in the Jahn- Teller distortion and the pairs 
of atoms forming the two dimers get closer to each other 
and form new weak Ga-Ga bonds with a length equal 
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to the intradimer distance. Finally a tetramer is formed 
where the four Ga atoms are equidistant and fivefold co- 
ordinated (three covalent Ga-As bonds and two weaker 
Ga-Ga bonds). The tetrahedron formed by the vacancy's 
first neighbor has 4 short bonds and 2 long bond, as can 
be clearly seen from Table UTTl 

This type of relaxation, "resonant bond" model, was 
first seen in calculations for the singly negative divacancy 
in silicouf^SiSS then for divacancies in GaAsi^i More re- 
cently, this relaxation pattern have been observed exper- 
imentally then confirmed using ab-initio cluster calcula- 
tions^ for the As Si — Vsi pair in silicon. 



2. Energetics 

We calculate the formation energies for all possible 
charge states of the As vacancy; the most relevant de- 
fects are reported in the first column of Table IIIII Most 
of the earlier calculations deal with the positively charged 
vacancy and do not go beyond the — 1 charge state; there 
are no recent calculations that report formation energies 
for the —2 and —3 charge states. 

In a 32-atom supercell LDA calculation of the forma- 
tion energy of charged defects, Northrup and Zhang find 
that the relaxed V^^ is stable for both intrinsic and p- 
type materials, ^•'^ with a formation energy 1.7 eV lower 
than the (-1-1) charge state. Our results are in complete 
disagreement with this calculation. The formation en- 
ergy we obtain does not favor the (-1-3) charge state un- 
der any doping or growing conditions: Ej:{V^^) =3.5 eV 
compared to Ej:{V^^)= 2.8 eV (in agreement with other 
calculations: £;^(V//)= 2.97 eV-^? using LDA and 3.09 
eVi^ using SCC-DFTB). This confirms that +1 has the 
lowest energy among all other defect. 

Figure 21 displays the formation energies as a function 
of Fermi level for various charge states of As vacancies at 
K. For p-type GaAs (Fermi level close to the valence 
band maximum) the -1-1 charge state has the lowest en- 
ergy, in agreement with recent results from Chadii^ 

We note, moreover, that lines for the neutral and the 
negative charge state intersect before those for the pos- 
itive and neutral state; the (0/—) ionization level is lo- 
cated well below the (-I-/0) level and represents a net 
signature of level inversion or of the so-called negative- 
U effect. It is therefore energetically more favorable to 
transfer two electrons at the same time to the defect level 
from the Fermi level with the reaction + 2e ^ V^^ . 

Such transfer is associated the strong Jahn- Teller dis- 
tortion discussed in Section IIV B II as the systems goes 
from the -t-1 to the —1 state. The negatively charged va- 
cancy remains stable for intrinsic and n-type GaAs but is 
superseded by the triply negative vacancy in the heavily 
n-doped GaAs corresponding to Ep = 1.27 eV (level b in 
figure^l): the (— /— 2) and (—2/— 3) defect levels almost 
collapse, thus the direct (— /— 3) transition is favored. 
Interestingly, the transition (— /— 3) is associated with a 



structural change from the pairing mode relaxation to the 
resonant mode relaxation as discussed in Section HV B II 
The four added electron are paired two by two as the 
energy gained from the structural relaxation overcomes 
the Coulombic repulsion for each of the two electrons, 
supplying a net effective attractive interaction between 
the electrons. 

Only the -1-1, and —1 charge states seem to have 
been studied previously by ab initio calculations. While 
Northrup et Zhang2& predict that arsenic vacancy in 
GaAs exist in the + charge state only, other calculated 
ionization levels split into two main categories classified 
in Table IVIl (i) A direct transition from a charge state 
to another is possible with only one electron transferred 
at time, and (ii) a negative-U effect for the (-1-/—) tran- 
sition. The first type of transition is found by Seong and 
LewiS)2^ who predict Vas to exist only in the + charge 
state, at 1.41 eV above the valence band this charge state 
changes to neutral then to negative charge state at the 
limit of the conduction band maximum. In an earlier 
study for unrelaxed As vacancies reported by Jansen and 
Sankey^^ ionization energies are located in the range of 
experimental band gap and located near the conduction 
band maximum. A negative-U effect is reported by both 
Cheong and Changi^ and Poykko et aS^ for the {+ / — ) 
transition, the transfer of the two electron occurs above 
the middle of the band gap. This negative-U behavior 
is confirmed by a recent calculation from ChadiSi which 
finds that the direct transition + 2e —f V^^ is favored 
after a Jahn- Teller distortion and that —1 charge state is 
the most stable for Fermi levels above midgap. 

Our results agree partially with this results. There is 
a negative-U effect but the level for the (-1-/—) transition 
is shallower. This might be due to the important struc- 
tural relaxation that affects the neutral and — charged 
vacancies. Moreover, as discussed above, we also find a 
second negative-U transition level at Ep = 1.27 eV with 
the reaction V^^ +2e ^ V^^ . 

Real GaAs crystals are far from being perfectly sto- 
chiometric, during growth there will be excess of Ga or 
As ions. A more general study concerns the effect of the 
growing conditions on the stability of the defect under 
certain doping conditions. In Figure [S] the various pan- 
els show the progressive doping of the GaAs sample and 
how the stability of the defect gets affected consequently. 
The Fermi level ranges between the valence band maxi- 
mum and the midgap where most of the ionization levels 
computed in this work have been identified. Each of the 
panels shows the transition between two charge states for 
As and Ga vacancies at the critical values of the ioniza- 
tion levels. For example at Ef=0.05 eV the first transi- 
tion for gallium vacancies take place, the two line collapse 
and are indicated as VqJ~^\ 

For p-type GaAs, at the As-rich limit in GaAs, corre- 
sponding to An/AHf close to —1, the dominant charge 
state are the As vacancies that probably compete with As 
antisites, while for the Ga-rich limit Ga vacancies have 
lower formation energies. For a Fermi energy at midgap 
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P-type (Ep =0 eV) Ep =0.4 eV 




FIG. 5: (Color online) Formation energies of Ga (solid line) and As (dashed line) vacancies in GaAs as a function of the growth 
conditions (A/x). Different panels are for different critical values of the Fermi level or ionization levels identified earlier. 



and for n-type GaAs, regardless of the growing condition, 
the triply negative charge state is the most stable among 
others and have the lowest formation energy. 

V. CONCLUSIONS 

We have presented a complete description of the en- 
ergetic and the relaxation geometry for relevant vacancy 
type defects in GaAs using the SIESTA ab initio pro- 
gram. Various convergence tests show that size effects, 
the completeness of the basis sets and the sampling of 
the Brillouin zone can become very important when the 
symmetry of the defect is broken or when the defect is 
highly charged. 

Using the DZP basis set, with a 216-atom unit cell and 
a 2 X 2 X 2 fc-point sampling, we find that Ga vacancies 



have shallow ionization levels below midgap in agreement 
with experiment, and do not show any recombination 
of the dangling bonds as was shown in earlier calcula- 
^jQjji0j25_ Pqj. -j-j^g j^ggg g^mjied As vacancy, we find that 

the ionization level (+/— ) of As is located in the lower 
half of the band gap and lies near the valence band while 
the second negative-U level (— / — 3) is located above 
midgap with a significant difference in the relaxation pat- 
tern reported earlier. The triply negative charge state for 
As vacancy reconstruct in the resonant bond mode, in a 
similar fashion as divacancies in Si and GaAs reported 
earlier. Finally, we find that only a few vacancy types 
can act as vehicles for self-diffusion of dopants in real 
GaAs devices under different doping and growing condi- 
tions, suvy as the triply negative Ga vacancy for intrinsic 
and n-type GaAs. 

These results will be used as a starting point for a 
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detailed study of self-difFusion using the SIEST-A-RT 
method presented elsewhere^ 
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